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Abstract. We present an analytic formalism describing structural properties of random uncorre- 
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between results of numerical simulations and our analytical predictions. 
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INTRODUCTION 

During the last years, there has been noticed a significant interest in the field of complex 
networks and a lot of interdisciplinary initiatives have been taken aiming at investiga- 
tions of these systems [1, 2, 3, 4]. It was observed that despite network diversity, most of 
real web-like systems share three prominent characteristics: small- world property (i.e. 
small average path length), high clustering and scale-free degree distribution [5, 6]. The 
observed universal properties let one understand that networks are not simply a sum 
of nodes connected by links. Nowadays, research on networks points mainly to the so- 
called emergent properties, that is systems global features and capabilities which are not 
specified by network design and are difficult or impossible to predict from knowledge 
of its constituents. As a result, the interest in topological characterization of real net- 
works gives way to growing interest in dynamical processes defined on such systems 
[7, 8]. We have already known how networks grow and how that growth process influ- 
ences network topology i.e. pattern of connections [9, 10]. We have also gained a certain 
understanding of how network structure mediates different transport phenomena like: 
disease (rumor) transmission in social networks [11] or information processing (virus 
infection) in computer networks [12]. However, there is still a number of open question 
like: what are the most efficient and robust topologies for different web-like systems 
[13]? How to effectively fight against criminal networks [14]? How to eliminate traffic 
jams in the Internet and other transportation systems [15]? 

This paper is devoted to uncorrelated random networks [16] (also known as random 
graphs or configuration model) which have been repeatedly shown to be very useful in 
modeling different phenomena taking place on networks. Although a number of other 



network models have been proposed, it was also proved that in many cases there are no 
qualitative differences between phenomena defined on both random graphs and those 
more sophisticated models. For example, the main results on percolation transition 
[17, 18] and epidemic spreading [12, 19] that were obtained for random uncorrelated 
networks are still valid for correlated and even non-equilibrium systems. The above 
certifies that despite its simplicity, configuration model is not merely mathematical toy. 
For the mentioned reasons, the issue of structural properties of random graphs and their 
reliable, mathematical description seems to be very important. It does not only excel in 
academic challenges but it can be helpful in further exploration of networks, including 
dynamical processes defined on discrete, web-like topologies. 

Following the idea, in this paper we proceed with microscopic description of the con- 
sidered networks and ask the most basic question related to their structural properties 
i.e. what is the probability that the shortest distance between two given nodes i and 
j is not longer than x. Basing on simple, intuitive arguments we derive the probabil- 
ity and apply it to calculate the main network characteristics like: the position of the 
phase transition at which a giant component first forms [20, 21], the mean component 
size below the phase transition [16], the size of the giant component [22, 17] and the 
average path length above the phase transition [23]. Although most of the enumerated 
macroscopic properties were previously calculated by means of other methods, espe- 
cially thanks to generating function formalism developed by Newman [16], we think 
that our approach is conceptually much simpler. Furthermore, since our derivations take 
advantage of primary concepts they do not represent a closed project but can be thought 
of as a background for further network investigations, particularly those associated with 
network dynamics. 

RANDOM UNCORRELATED NETWORKS - MICROSCOPIC 

DESCRIPTION 

Random uncorrelated network with arbitrary degree distribution P(k) is the simplest 
network model. In such a network the total number of vertices /V is fixed. Degrees of all 
vertices are independent, random integers drawn from a specified distribution P(k). One 
can easily create the network following the below instructions: 

i. Prepare N nodes i = 1,2, . . . ,N. 

ii. Attach to each node k[ ends of connections taken from the given distribution P(k). 

iii. Connect at random ends of connections. 

The above procedure provides maximally random, uncorrelated networks. From the 
mathematical point of view the lack of correlations means that the probability P(kj/kj) 
that an edge departing from a vertex j of degree kj arrives at a vertex i of degree ku 
is independent of the initial vertex j. The above translates into the fact that the nearest 
neighborhood of each node is the same (in statistical terms). Provided that the degree 
distribution in uncorrelated network is P(k) one can show that the degree distribution of 
nearest neighbors, that is equivalent to the probability that an arbitrary link leads to a 



node i of degree kj, is given by 



P(k l /kj)=Q(k l ) = ^ ) P(k i ), (1) 
and respectively the average nearest neighbor degree equals 

= (2) 

At the moment, before we proceed with our derivations, let us introduce notation that 
will be used throughout the paper: 

i. p\- (x) denotes the probability that there exists at least one walk of length x between 
two given nodes i and j. If such a walk exists then the shortest path between the 
two nodes is < x. 

ii. pjj(x) gives the probability that none among walks of length x occurs between the 
two nodes, and respectively 

p7 j { X ) = \-p±{x). (3) 

iii. p*j(x) describes the probability that the shortest distance between i and j is equal 
to x. In the limit of dense networks (above percolation threshold) one can deduce 
that if there exists at least one walk of length x between i and j then walks longer 
than x also exist. It follows, one can assume that pfj{x) (see i.) expresses also the 
probability that the shortest distance between the two nodes is not longer than x and 

p* j (x)=p±(x)-p±(x-l). (4) 

Now, we try to derive formulas describing the above probabilities. At the beginning, 
let us imagine a random walker that starts at a given node i and assign the walker one- 
step memory in order to avoid backward steps. Next, let the walker perform x steps. It is 
easy to see that in average the walker may pass 

w(i,x) = k i ((k nn )-iy- 1 (5) 

different walks. If, at least one such a walk ends at the node j one can say that both 
vertices i and j are at most x—ih. neighbors i.e. the shortest path between the two nodes 
is not longer than x. It follows that 

/W(i,x) 

p!M)=P\ U w£j]> ( 6 ) 



where W-i+j represents a single walk between i and j, whereas the sum goes over all 
possible x— walks starting at i (5). 



To simplify the last expression (6) let us assume mutual independence of all walks 

(r) 

Wj 'j. Then, the considered formula may be rewritten in the following form (see Ap- 

*• J 

pendix A) 



Ptj( x ) -1-exp 



W(i.x) 



(7) 



To finalize derivation of pfj{x) note that, due to the lack of correlations, every walk 
starting at i has the same chance to end up in j. The probability of such an event is just 
equal to the probability that the last connection in the sequence of edges representing 
the considered walk is attached to j 



P W; 



ft) 



(k)N 



(8) 



Inserting the last formula into (7) and taking advantage of both (5) and (2) one gets 

p+(x) = l-F(x), (9) 

where 

r kjkj f(k 2 ) ^ 
(k)N V (k) 



F(x) = exp 



x-l 



(10) 



The last formula for pfj(x) constitutes the most important result of the paper. Let us 
stress that it does not include any free parameters, therefore one can directly employ it to 
calculate all the interesting structural characteristics of the considered networks. In the 
following sections we show, how the formula for percolation threshold in random graphs 
several times naturally emerges from our approach. We also calculate, in turn, the mean 
component size below the percolation transition, the size of the giant component and 
also the average path length above this transition. 



MEAN COMPONENT SIZE BELOW PERCOLATION 

THRESHOLD 

To calculate the mean component size below the percolation threshold one can make use 
of the introduced probabilities pfj (x) and pj- (x) . 

Taking advantage of (9), the probability that none among the walks of length x 
between the nodes i and j exists is given by (3) 

pr.(x) = exp 

and respectively the probability that there is no walk of any length between these vertices 
may be written as 



:)N { (k) J 




The value of P- strongly depends on the common ratio of the geometric series present 
in the last equation. When the common ratio is greater then 1, i.e. (k 2 ) >2(k), random 



graphs are above the percolation threshold. The sum of the geometric series in (12) tends 
to infinity and Pfj = 0. Below the phase transition, when (k 2 ) <2(k), the probability that 
the nodes i and j belong to separate clusters is given by 



p ij = ex P 



kikj 



1 



(13) 



N (2(k)-(k 2 ))_ 

and respectively the probability that i and j belong to the same cluster may be written as 



P+ 

ij 



1-Pij = 1-exp 



kikj 



1 



N (2(k)-(k 2 )) 



(14) 



Now, it is simple to calculate the mean size of the cluster that the node i belongs to. It 
is given by 



(s)(k i ) = i+i,p(k j )p+~ 



(k) 



2(h) - (k 2 ) 



ki + l. 



(15) 



Note, that the mean size of the component that a node i belongs to, is a linear function 
of degree k[ of the node (see Fig. 1). The last transformation in (15) was obtained by 
taking only the first two terms of power series expansion of the exponential function in 
(14). Averaging the above expression (15) over all nodes in the network one obtains the 
known formula [16] for the mean component size in random graphs below the phase 
transition 



2(k) - (k 2 ) ■ 



(16) 



As in percolation theory [24], the mean cluster size diverges at 

(k 2 ) = 2(k), (17) 

one more time certifying that the expression (17) describes the position of the perco- 
lation threshold in random uncorrelated networks with arbitrary degree distributions 
[20,21]. 



SIZE OF THE GIANT COMPONENT 

When (k 2 ) > 2{k) the giant component (GC) is present in the graph. Its relative size S, 
i.e. the probability that a node belongs to GC, is an important quantity and is often iden- 
tified as the order parameter of the percolation transition. In this section, we demonstrate 
how to calculate the size of the giant component S in uncorrelated networks with arbi- 
trary degree distributions. The underlying concept of our derivations is closely related 
to the method of calculating S in Cay ley tree and originates from Flory (1941) [25, 24]. 

At the beginning, we deal with classical random graphs of Erdos and Renyi, then 
we generalize our calculations for the case of random graphs with arbitrary degree 
distributions and we show that our derivations are consistent with the formalism based 
on generating functions that was introduced by Newman et al. [16]. 



Classical random graphs of Erdos and Renyi (ER) 



In general, one defines the classical random graph as N labeled nodes and every pair 
of the nodes being connected with probability p [26] . The probability that an arbitrary 
node i belongs to the giant component is equivalent to the probability that at least one 
of its /V — 1 possible links connects it to GC. In order to describe the above equivalence 
by means of mathematical expression, let us assume that A^jy represents an event: the 
connection {i,j} (if exists!) leads to the giant component. The notation A^-jj. refers to 
each of N(N — l)/2 possible links, not only to those existing! Now, the size of the giant 
component is given by 

and due to the mutual independence of different links, the last formula may be rewritten 
as (see Appendix A) 

r n-\ 

5=1 -exp • (19) 

Next, to further simplify the expression for S note that the mentioned mutual indepen- 
dence implies that every link has the same probability to belong to the infinite cluster 

V{iJ} P (A {iJ} )=R. (20) 




<k> 



FIGURE 2. Size of the giant component S versus the average degree (k) in classical random graphs of 
size N = 10000. The scatter plot represents numerical data whereas the solid line gives the solution of the 
Eq.(23) 



Inserting (20) into (19) one gets 

S=l-exp[-(N-l)R]. (21) 

To finalize derivation of S one has to find R. To do so, let us recall that R describes 
the probability that an arbitrary node i is connected to the giant component through a 
fixed link {/, j}, where j is another arbitrary node. Now, if j belongs to GC it means that 
at least one of its connections also leads to the giant component. Given that every node 
may have N—l links the formula for R may be expressed as a product of the probability 
of a link {i,j} and the probability that at least one of N — 1 possible links emanating 
from j connects the considered node to the giant component. Summarizing the above 
considerations one obtains a self-consistency equation for R 

R = nA {hj] )=pp([}A {hm A (22) 

\m=l / 

= p(l-exp[-(N-l)R}). 

Finally, comparing both relations (21) and (22) it is easy to see that R = pS and the 
expression for the giant component in classical random graphs may be rewritten in the 
following form [22] 

S=l-exp[-(*)S], (23) 

where (k) = pN (see Fig. 2). At the moment let us point a certain interesting property of 
the last equation, that makes the equation very intuitive example of percolation transition 
(at least for those acquainted with Ising model). Below the percolation threshold (i.e. for 
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FIGURE 3. Graphical solution of the Eq. (23). 



(k) < 1 [ ) the identity has only one solution 5 = (see Fig. 3). Above the threshold 
another solution 5^0 appears signifying transition of the system to the ordered phase. 

Giant component in random graphs with arbitrary degree 

distributions 

In this section, taking advantage of the intuition gained from the analysis of the giant 
component in classical random graphs, we develop a more general approach allowing 
to calculate the size of the giant cluster in random networks with arbitrary degree 
distributions P(k). 

In the case of classical random graphs all nodes have been considered equivalent. 
It is not acceptable in the case of random graphs with a given degree sequence P(k), 
where every node has a fixed number of connections. In order to meet the requirements 
imposed by P(k), we have to slightly rearrange the previous meaning of the probability 
R (20). Now, let us assume that R describes the probability that an arbitrary but existing 
(!) connection belongs to the giant component. It is also useful to think of R as the 
probability that following arbitrary direction of a randomly chosen edge one arrives at 
the giant component. In fact, we know that following an arbitrary edge we arrive at a 
vertex of degree k. The probability that the node, we have just arrived at, is connected to 
GC equals 1 — ( 1 — R) k ~ l . The last relation simply expresses the probability that at least 
one of k — 1 edges emanating from i and other than the edge we arrived along connects i 
to the giant component 2 . Summarizing the above considerations makes simple to write 



Degree distribution in classical random graphs is Poissonian, therefore the expression for percolation 
threshold (17) simplifies to (k) = 1. 

2 Here, we do not take advantage of the Lemma 1 because it works well only in the limit of large number 
of the contributing events n » 1 . In the case of small n the error q of the Lemma 1 can not be neglected. 



the self-consistency condition for R (compare with Eq. (22)) 



k 



(24) 



where Q(k) is given by (1). Then, knowing R it is easy to calculate the relative size S 
of the giant component that is equivalent to the probability that at least one of k links 
attached to an arbitrary node connects the node to GC 



To make derivations of this section more concrete, we should immediately introduce 
some examples of specific networks. Since however, one can show that both formulas 
(25) and (24) are completely equivalent to equations derived by other authors (see 
Appendix B), we just refer the reader to analyze examples presented in those related 
papers [16, 27]. 



AVERAGE PATH LENGTH IN RANDOM UNCORRELATED 

NETWORKS 



We turn now to the quantitative description of the small-world effect in random un- 
correlated networks with arbitrary degree distributions P(k). To our knowledge, the be- 
low derivations are the simplest and the most accurate among those previously reported 
[16, 28, 23], therefore we illustrate them with a larger number of examples. 

Taking advantage of (9), we are in a position to calculate p*j(x) (4) i.e. probability 
distribution for the shortest distance between any two nodes i and j 



where F(x) is given by (10). Although, due to the condition p*j(x) > the last expression 
is only correct above the percolation threshold, the formula (26) is very important and, 
given that the giant component contains almost all vertices, it allows to find a number 
of structural properties of the considered networks. For example, averaging (26) over all 
pairs of nodes one obtains the intervertex distance distribution p(x) = ({p*j(x))i) j. It is 
also possible to calculate z x i-e. the mean number of vertices a certain distance x away 
from a randomly chosen node i. The quantity can be obtained as z x = N J p* i j{x)P{kj)dkj. 
At the moment, let us note that taking only the first two terms of power series expansion 
of both exponential functions F in (26) one gets the relationship z x = Zi(z2/z\) x ~ l = 
{k)({k nn ) — l) x_1 (compare it with (5)) that was obtained by other authors given the 
assumption of a tree-like structure of random graphs. Finally, the average path length 
(APL) between any two nodes i and j of degrees k[ and kj is simply the expectation 
value of the distribution (26) 




(25) 



pUx) = F(x-l)-F(x), 



(26) 



oo 



hj(ki,kj) = Y, x Pij( x ) = E F W- 



(27) 
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FIGURE 4. Average path length Ier versus network size N in ER classical random graphs with 
(k) = pN = 4, 10,20. The solid curves represent numerical prediction of Eq. (31). 



The Poisson summation formula allows one to simplify the above sum (see Appendix C) 

t . * -ln(^)+ln((fc 2 )-(fc))+lnjV-y 1 

lij{khkj) ~ ln«*2)/<*)-l) + 2' (28) 

where f— 0.5772 is the Euler's constant. The average intervertex distance for the whole 
network depends on a specified degree distribution P(k) 

-2(\nk)+\n((k 2 )-(k))+\nN-r 1 

ln«*2>/(*>-l) + 2- ^ J 

As one could expect, the both formulas (28) and (29) diverge when the considered 
networks approach percolation threshold (k 2 ) = 2(k) (17). 

To test the formula (29) we have performed numerical simulations of the well-known 
networks: classical random graphs proposed by Erdos and Renyi (ER), single-scale 
networks with exponential degree distribution (EXP) and scale-free networks (SF). 

Classical random graphs (ER). For these networks the degree distribution is Pois so- 
man 

P(k) = ~ 1 ^ L - (30) 

However, since (In A:) cannot be calculated analytically for Poisson distribution thus 
the APL may not be directly obtained from (29). To overcome this problem we take 
advantage of the mean field approximation and assume that all vertices within a graph 
possess the same degree Vjki = (k). It implies that the APL between two arbitrary nodes 
/ and j (28) also describes the average intervertex distance characterizing the whole 
network 

lniV-y 1 

'-=M^) + 2- (31) 
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FIGURE 5. Characteristic path length Iexp versus network size N in exponential networks. Solid lines 
represent Eq. (33). 



Until now only a rough estimation of the quantity has been known. One has expected 
that the average shortest path length of the whole ER graph scales with the number 
of nodes in the same way as the network diameter. We remind that the diameter d of 
a graph is defined as the maximal distance between any pair of vertices and <1er = 
lniV/ln(piV) [3]. Figure 4 presents the prediction of the equation (31) in comparison 
with the numerically calculated APL in classical random graphs. 

Exponential networks (EXP). Now, let us suppose that the degree distribution is 
exponential 

e (l-k)/a 

P(k) = , (32) 

a 

where k = 1,2,... 3 . Applying the distribution to Eq. (29) immediately provides the 
formula for the average path length in the considered single- scale networks 

\nN + ln(2a 2 + a)-2e 1 / c T(0,l/a)-y 1 

1 exp = wo 2 ! — ^ — n — rr\ ( 33 ) 

ln(2a / + a) — m(a + 1) 2 

where T(0, 1 /a) is an incomplete gamma function. Figure 5 shows that the obtained 
formula perfectly fits numerical results obtained for different values of the parameter a. 

Scale-free networks (SF). As mention at the beginning of the paper degree distribu- 
tions are scale-free in most of real systems 

(a- 1 W" -1 ) 
p(k) = (a !)m m 



3 We start to enumerate degrees from k = 1 instead of k = in order to prevent construction of very sparse 
networks with a large number of separated nodes. 



where k = m,m+ \,...,k max . It was found that the most generic mechanism driving 
real networks into scale-free structures is the linear preferential attachment. The sim- 
plest model that incorporates the rule of preferential attachment was introduced by 
Barabasi and Albert [6]. Other interesting mechanisms leading to scale-free networks 
were proposed by Goh et al. [29] and Caldarelli et al. [30]. Unfortunately, the mentioned 
mechanisms leading to scale-free network topologies incorporate structural correlations, 
making our approach useless. The simple procedure of generating random uncorrelated 
networks that was described at the beginning of this paper also fails when applied to 
fat-tailed degree distributions with diverging second moment (k 2 ) . In particular, the pro- 
cedure may not be applied to generate uncorrelated scale-free networks (34) with the 
scaling exponent 2 < a < 3. In order to avoid the inconvenience connected with those 
scale-free distributions we apply the so-called structural cutoffs by imposing the largest 
degree to scale as k max = a/ (k)N [31, 32, 33], instead of unbounded cutoff predicted by 
extreme value theory k max ~ AT 1 /( a_1 ). 

We have found that depending on the value of the exponent a one can distinguish 
three scaling regions for the average path length (29) in scale-free networks (see Fig. 6). 
In the limit of large systems N — > °°, APL scales according to relations 4 

• for a > 3 

Z«> 3 ~ lnJV, (35) 

• for a = 3 



• for 2 < a < 3 

Note that although the results for a > 3 are consistent with estimations obtained 
by other authors [23, 34], the case of 2 < a < 3 is different. In opposite to previous 
estimations signaling the double logarithmic dependence l a<3 ~ \n\nN, our calculations 
for the same range of a predict that there is a saturation effect for the mean path length in 
large networks. Unfortunately, at the moment it is impossible to assess what is the correct 
estimation (to perform reliable tests very large networks, even with N > 10 10 , should be 
analyzed). On the other hand, it is truth that looking at Figure 6 one can observe two 
effects: i. the difference between results of numerical calculations and our analytical 
predictions continuously grows with N, ii. denser networks are better described by our 
approach. The first effect may result form the fact that our mean-field derivations work 
worse for networks with degree distributions characterized by large fluctuations (note 
that the second moment of scale-free distribution described by the exponent a < 3 
increases with AT), whereas the second effect may be related to the fact that our approach 
(in particular Eq. (26)) works better for networks further above the percolation threshold 



4 The exact formulas for APL in those regions are very complicated therefore we have decided to omit 
them. 



Z«= 3 ~ ™- (36) 



l a < 3 = const. (37) 
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FIGURE 6. Average path length versus number of nodes N in scale-free networks with a =4 (CI), 
a = 3 (C2) and a = 2.5 (C3). The scatter data represent numerical calculations. Solid curves with open 
symbols express analytical predictions of Eq. (31). 



CONCLUSIONS 



To conclude, in this paper we have presented theoretical approach for metric properties 
of uncorrected random networks with arbitrary degree distributions. We have derived a 
formula for probability pfj(x) (9) that there exists at least one walk of length x between 

two arbitrary nodes i and j. We have shown that given pfj(x) one can find a number 
of structural characteristics of the studied networks. In particular, we have applied our 
approach to calculate the mean component size below the percolation transition, the 
size of the giant component and the average path length above the phase transition. We 
have successfully applied our formalism to such diverse networks like classical random 
graphs of Erdos and Renyi, single-scale networks with exponential degree distributions 
and uncorrelated scale-free networks with structural cut-offs. In all the studied cases 
we noticed a very good agreement between our theoretical predictions and results of 
numerical investigation. 



APPENDICES 
Appendix A 

Lemma 1 If A i , A2, . . . , A n are mutually independent events and their probabilities fulfill 
relations V/P(A,-) < £ then 

P(|jA ! ) = l-exp(-f P(Ai))-Q, (38) 
(=1 (=i 

where 0<Q< U]H(n£)j / j\ - (1 + e) n . 

Proof. Using the method of inclusion and exclusion [35] we get 

^(U A <) = L(-l) ;+1 S(j), (39) 
i=i 7=1 

with 

S(j)= £ P(A h )P(A i2 )...P(A ij ) 

\<i\<i2<...<ij<n 

£p{A t )\ -Qj, (40) 



1 



J'! 



where < Qj < (n ] / j\ — £ J - The term in bracket represents the total number of 

redundant components occurring in the last line of (40). Neglecting Qj it is easy to see 
that (1 — P(UAi)) corresponds to the first (n+ 1) terms in the MacLaurin expansion 
of exp( — £P(A,-)). The effect of higher-order terms in this expansion is smaller than 



R < (ne) n+l I (n + 1)!. It follows that the total error of (38) may be estimated as Q < 
L"=i Qj +R- This completes the proof. 

Let us notice that the terms Qj in (40) disappear when one approximates multiple 
sums Li<i 1 <i 2 <...</-< w by corresponding multiple integrals. For £ =A/n <C 1 the error 
of the above assessment is less then A 2 exp(A)/n and may be dropped in the limit n — > °°. 



Appendix B 

It is easy to show that our formulas (25) and (24) are completely equivalent to 
equations derived by Newman et al. [16] by means of generating functions 

S=l-G (v), (41) 

where v is the solution of equation given below 

v = Gi(v). (42) 

We recall that Gq{x) is the generating function for the degree distribution 

Go(*) = !>(*)**, (43) 

k 

whereas G\(x) is related to Q(k) (1) 

Let us start with Eq. (24) that may be transformed in the following way 

R = £ Q(k) - £ Q(k) ( 1 - R) k - 1 = 1 - Gi ( 1 - 7?) . (45) 

yt it 

Note, that the last formula directly corresponds to Eq. (42) with 

v=l-R. (46) 

Taking into account the last identity one can show that the expression (25) may be 
transformed into Eq. (41) in a similar way. Now, it is clear that the unknown parameter 
v in both Eqs. (41) and (42) has the following meaning - it describes the probability that 
an arbitrary edge in a random graph does not belong to the giant component. 



Appendix C 

The Poisson summation formula states 

oo 1 

£F(x) = -F(0)+ (47) 

x=0 Z 



POO L * J / POO 

/ F(x)dx + 2'Y[ I F(x)cos(2n7lx)dx 
Jo „ti V> 



Applying the formula to (28) 



lij(ki,kj) = £exp 
one realizes that in most of cases 



m-(k))N\ (k) 



kikj 



k(k j 



m-(k))N 

that gives F(0) = 1. One can also find that 



~0 



Jo 



F{x)dx = —Ei 



kjk j 



((k 2 )-(k))N 



/In 



-1 



(48) 



(49) 



(50) 



where Ei(y) is the exponential integral function that for negative arguments is given by 
Ei{— y) = y + lnj [36], where f — 0.5772 is the Euler's constant. Finally, it is easy to 
see that owing to the generalized mean value theorem every integral in the last term of 
the summation formula (47) is equal to zero. It follows that the equation for kj{k^ kj) is 
given by (28). 
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